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INTRODUCTION 

This is the final report on Atmospheric Sciences Section, National 
Science Foundation Grant GA~524, a two-year research project designed to 
investigate the winter thermal budget of the Great Lakes (principally Lake 
Michigan), as a necessary sequel to a previous project supported by the 
Atmospheric Sciences Section of the National Science Foundation under NSF 
Grant GP-2411, titled "Growth of Ice on Lake Michigan." The final report 
of the previous project was published in 1966 as Special Report No. 26 of 
the Great Lakes Research Division of the University of Michigan (Heap and 
Noble 1966) o The first annual report of this project was published as 
Special Report No. 32 of the Great Lakes Research Division (Noble and Ewing 
1967) o A related study of the summer temperature structure of Lake Michigan 
for 1963-66 was carried out under the cognizance of the current project, 
and was reported in Special Report No. 30 of the Great Lakes Research Divi- 
sion (Noble 1967). 

Inasmuch as this is the final report for the research project, the 
previous results will be summarized, as necessary, for completeness. New 
data and analyses will be presented in detail. Readers may obtain copies 
of previous reports from the Great Lakes Research Division for full data 
reports published earlier. 

Two graduate students have participated heavily in this project and 
have made significant contributions to the results of the program: Joseph C. 
Huang, who has done much of the computer programming and who is doing a 
related dissertation problem, and Karen J. Ewing, co-author of this report. 
Of the three undergraduate students working on the project, April Michaelis 
has made a most valuable contribution and has co-authored a related data 
report. Special Report No. 39 of the Great Lakes Research Division (Noble 
and Michaelis 1968). 



ICE PREDICTION INDICES 

The research program supported by this grant was designed to verify 
potential ice prediction techniques indicated by a previous study reported 
by Heap and Noble (1966) o 

Within constraints to be described in following sections of this report, 
and within verification conditions imposed by the somewhat limited ice cover 
observations for Lake Michigan, an operational ice prediction technique is 
promulgated for initial implementation by the UeS. Coast Guard, U^S* Weather 
Bureau, and the Lake Carriers Association « The operational method is ini- 
tially designed for Lake Michigan, with the understanding that the basic 
methods may be subsequently applied to all of the Great Lakes as substantia- 
ting data are developedo It is suggested that the technique be extended to 
the remainder of the Great Lakes with the cognizance and cooperation of the 
Canadian Department of Transport o 

loe 'Potential Index 

As indicated in the previous study, and as verified for a five-year 
period (Noble and Ewing 1967), the early-winter temperature structure of 
Lake Michigan may be used as an indicator of the ice-forming potential of 
Lake Michigan o 

If the temperature of the mid-lake water mass (the deep portions of 
the southern basin of the lake) is of the order of 9*^C for the first week 
of December^ there is sufficient heat storage in the lake to reduce ice 
formation to a nearly insignificant amount (eogo the winter season of 1963- 
64) o If the temperature of the mid-lake water mass is of the order of 6°C 
for the first week of December, ice formation may approach 100% coverage of 
the lake basin (eogo the winter season of 1962-63). These criteria were 
established from weekly averages calculated from U.S. Public Health Service 
Federal Water Pollution Control Administration buoy stations in Lake Michigan. 

BT data from the USCGC WOODBINE, taken during her late November-early 
December run to recover aids to navigation in the southern basin of Lake 
Michigan^ have served as a useful approximation to average thermograph data 
for an Ice Potential Indexo Further ^ since the water column of the lake is 
nearly isothermal at this time of year, it has been demonstrated that accur- 
ately calibrated surface temperature measurements taken from the WOODBINE 
may be used as the Ice Prediction IndeXo (An accurate, calibrated, ther- 
mistor bridge was installed in the main sea chest of the WOODBINE in the 
fall of 1967 o) Table 1 summarizes the Ice Potential Index data for the 
winter seasons 1962-63 through 1967-68 o 

The following technique is suggested for the determination of an annual 
Ice Prediction Index for Lake Michigan^ 

A calibration check should be run on the WOODBINE'S thermistor bridge 
during early November o During the fall cruise to recover aids to navigation 
(typically Grand Haven - Milwaukee - Chicago - Grand Haven during the last 
week of November and the first week of December) , the WOODBINE should read 



Table lo Ice Potential Index data for winter seasons 1962-63 through 1967-68, 

1962-63 Average water temperature 5o9°C, 2 December-9 December 1962, 
buoy station 18 (43°08'N, 57°24o5'W)« 
Ice potential: Heavy o Ice cover: Record ice formation. 

1963-64 Average water temperature 8o8°C, 2 December-9 December 1963, 
buoy station 18 (43°08'N, 57''24o5'W)c 
Ice potential: Light o Ice cover: Very light ice formation, 

1964-65 No water temperature measurements « 

1965-66 WOODBINE BT's water temperature 7«3'*C to 9,4*'C, Southern Basin 
Lake Michigan, 3 December-9 December 1965 o 
Ice potential: Light o Ice cover: Very light ice formation, 

1966-67 WOODBINE BT's water temperature 4.2°C to 8o2°C, 3 December-9 
December 1966, Southern Basin Lake Michigano 
Ice potential: Moderate p Ice cover: Moderate to light, 

1967-68 WOODBINE BT's and thermistor bridge 3.5°C to 5«3°C, 13 December- 
18 December 1967, Southern Basin Lake Michigano 
Ice potential: Heavy o Ice cover: Early ice melted mid-winter. 
Later ice held to late in seasono 



surface water temperatures from her thermistor bridge (to the nearest 0,1**C) 
at each point of departure upon entering and leaving harbor entrances, and 
at each hour while away from the docko Each day (say at 1000 hours, depend- 
ing upon traffic) the WOODBINE should transmit a message to 9th District 
Coast Guard Headquarters giving the surface water temperature readings for 
the preceding 24 hours o The message would be coded in the following format: 

LLLL 1111 ttt, LLLL 1111 ttt, oo« 

where LLLL is the latitude in whole degrees, minutes, and tenths 
of minutes (omitting the ten-degree digit 4) 

where 1111 is the longitude in whole degrees, minutes and tenths 
of minutes (omitting the ten-degree digit 8) 

where ttt is the water surface temperature to Ool^Co 

Upon receipt of the coded message at Coast Guard headquarters in 
Cleveland, the WOODBINE temperatures would be given a subjective analysis 
to derive the Ice Prediction Indexo Subjective analysis is required to in- 
terpret the temperatures obtained around the southern basin of Lake Michigan 
in terms of a temperature representative of the deep-water mass of the lake. 
In other words ^ ship injection temperatures are used to approximate the 
weekly average, deep-water thermograph data obtained in 1962-63 and 1963-64. 
Equivalent deep-water temperatures (representative of the southern basin of 
the lake), between 6°C and 9°C for the first week in December, indicate a 
potential for the formation of an ice cover ranging from extremely heavy to 
very light, respectively o 



It is strongly suggested that all of the Coast Guard cutters on the 
Great Lakes be provided with precision water intake temperature indicators 
similar to that installed on the WOODBINE <^ and that similar statistics be 
developed for Ice Potential Indices for all of the Great Lakes «> 

Freeze-vjp Date and Rate, of lee Formation 

Verified statistics have not been developed ^ but a method similar to 
that of Richards (1964) has been indicated for the short-range (one week) 
prediction of the freeze-up date and the rate of ice formation of Lake 
Michigano This technique is based on the analysis of daily air temperature 
data from Milwaukee o As will be developed in a following section of this 
report J the method is to compute the mean of the high and low temperature 
for each day^ convert this mean to temperature expressed in degrees Celsius, 
and maintain a cumulative plot of the sum of the daily mean temperatures 
from 1 October of each year^ Since the freezing point is 0®C, the cumula- 
tive plot represents degree=-days of freezing (or thawing) exposureo The 
cumulative plots (Noble and Ewing 1967; Heap and Noble 1966) show a positive 
increase from 1 October until approximately 25 November ^ when a maximum is 
reached^ The curve then goes negative from the last week in November, in- 
dicating the onset of freezing weather o The date of the maximum of the 
curve is a precursor of the freeze-up date, and the rate of decrease of the 
cumulative curve indicates the rate of ice formationo The curve begins to 
decrease as negative degree-days (freezing exposure) are algebraically added 
to the cumulative curve o 

Inadequate ice observations over a sufficient period of record for Lake 
Michigan have prevented the establishment of empirical criteria for the 
definition of freeze-up dates or of the rate of ice formation. These data 
are currently being developed through systematic ice observations being 
carried out by the UoSo Lake Survey o 

The Ice Prediction Index provides an estimate of the amount of ice to 
be expected in any given seasono The cumulative Celsius degree-day plots 
provide an indication of the f reeze=up date, and the rate of formation of 
ice to the predicted maximumo 

The long-range warming or cooling trend of Lake Michigan may be esti- 
mated for four-year periods by the method of Ayers (1965) o 
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A PRELIMINARY STUDY OF THE LAKE MICHIGAN HEAT BUDGET 

Ayers (1965) suggested that annual averages of air temperatures, cloudy 
days, and storm passages over the Lake Michigan basin can be combined as 
moving 5-year means to predict warming and cooling trends of the water tem- 
peratures measured at the Chicago water intakes. These climatological cal- 
culations have been carried out from data available for the last century, 
and indicate the observed water temperature trends approximately five years 
in advance o Ayers' resultant predictor curve has been extended to 1965, 
and is reproduced in Figure 1. The predictor curve indicates that the lake 
should cool and remain at lower temperatures through approximately 1970 ^ 
Table 1 showed a cooling trend from 1965 through 1967. It is anticipated 
that the lake will remain at low temperatures for another two years. 

In order to study the heat storage in the lake basin, a computer pro- 
gram has been developed to analyze all BT data from the Great Lakes Research 
Division for monthly average water temperatures for 10-meter layers for each 
month for each one-degree latitude/longitude square in the Lake Michigan 
Basine The data have been published for 1963-1966 by Noble (1967). The BT 
punch card format used by the Great Lakes Research Division, the analysis 
program, and the data for 1954, 55, 61, 62, and 1967 (all available data) 
are presented in Appendix B. 

A second program has been developed that analyzes each individual BT 
cast, determines whether a thermocline exists, or whether the cast is iso- 
thermal, and defines the upper and lower limits of the boundary layer and 
the depth of maximum temperature gradient <» The analysis program and the 
summary data giving the limits of the boundary layer and the depth of maxi- 
mum gradient by month for each degree of latitude of Lake Michigan for 1954 
through 1967- are presented in Appendix C. 

Because of the lack of winter water temperature data (except for the 
1962-63 and 1963-64 seasons) and because of the lack of complete ice recon- 
naissance data for the whole Lake Michigan basin, the winter heat transfer 
from the water mass of Lake Michigan must be estimated from fall lake tem- 
peratures, spring lake temperatures, and the sensible heat transfer into 
the air mass as it crosses the lake under the effect of the prevailing 
westerly winds o This estimate of heat transfer does not provide adequate 
data for heat budget calculations, but may be used to estimate trends in 
the heat storage in the lake. In order to measure the solar radiation avail- 
able for lake heating during the winter, integrating solar radiometers were 
developed and tested in the Ann Arbor area during the winter of 1966-67, as 
described in the previous reports The integrating radiometers were modified 
slightly and placed at Coast Guard stations during the winter of 1967-68 at 
St^ Joseph, Muskegon, and Manistee, Michigan; Two Rivers, and Milwaukee, 
Wisconsin; and Chicago, Illinois, to determine the daily integrals of solar 
radiation incident upon the whole lake basin. Unfortunately, due to com- 
ponent quality control problems and to the inexperience of the operators at 
some stations, not enough solar radiation data were obtained to be useful 
for a thermal study a 

Daily air temperature records were examined for Milwaukee, Muskegon, 
Two Rivers, Manistee, Chicago, and Benton Harbor. The average between the 
daily high and low temperatures was calculated and converted to temperatures 
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expressed in degrees Celsius o Since the freezing point is 0*^C, negative 
values indicate freezing exposure, and positive values indicate thawing 
exposure o October 1 was arbitrarily defined as the beginning of the winter 
season^ and the daily average temperatures were cumulatively summed. Since 
the average daily air temperature was above freezing during the month of 
October, the cumulative degree-day curve increases to a maximum that occurs 
at the onset of freezing weather o The negative degree-days that are indi- 
cative of freezing weather are subtracted from the previous cumulative total. 
The cumulative curve reaches a minimum and starts to increase again with 
the onset of warm spring weather in April. The date of the maximum of the 
cumulative degree-day curve indicated the onset of freezing weather, and 
therefore, of the freeze-up dateo Figure 2 gives the cumulative degree-day 
curve for Milwaukee for the winter seasons of 1962-63 through 1967-68. The 
difference between the maximum and minimum values, of the cumulative degree- 
day curve is a measure of the severity (freezing exposure) of the winter 
seasono It is seen that the curve maximum occurs at dates varying from 25 
November to 16 December, and that the curve minima occur between 2 March 
and 31 Marcho Freezing weather begins at the time of the. curve maximum, 
and thawing weather at the time of the curve minimum. The dates of the 
maxima and minima of the curves are indicators of the freeze-up and break- 
up dates, respectively o Of course, accurate measurements of the solar radi- 
ation must be considered to obtain a precise prediction of the freeze-up and 
break«=up dates o 

Table 2 summarizes the cumulative degree-day data for three pairs of 
weather stations for the winter seasons of 1962-63 through 1963-64o The 
station pairs 3, Two Rivers - Manistee, Milwaukee - Muskegon, and Chicago - 
Benton Harbor were selected to be representative of the northern, middle, 
and southern portions of the lake basinp The dates and the cumulative 
degree-day values are given for each of the station records. The algebraic 
difference between the maximum and minimum values is the freezing exposure 
for each of the records o The Milwaukee records show that the severity of 
the winter seasons varied from 866 degree-days of freezing exposure in 1962- 
63 (the record ice year) to 439 during 1967-68. The very low freezing 
exposure during 1967-68 explains why, in spite of an ice potential index 
that indicated a possibility of extremely heavy ice cover, relatively little 
ice was formed as compared with the record year of 1962-63. On the other 
hand, the heavy ice potential index explains the persistence of the late- 
winter ice long into the spring of the year^ 

For the six stations considered, the winter freezing exposure ranged 
from 192 to 866 Celsius degree-days during the six winter seasons under 
study o The summer (7 April to 1 October) heating period showed approximately 
3000 Celsius degree-days for all station records « The data for freezing 
exposure and heating (thawing) exposure are to be considered only as exposure 
indices^, and cannot be construed as direct indicators of heat flux. For 
example, Table 3 gives the significant points for the cumulative degree-day 
curves from 6 April 1962 through 7 April 1968o The date points given are 
the beginning of the defined summer period, 7 April, the beginning of the 
defined winter period, 1 October, and the dates of the winter maxima and 
minima, approximately 10 December and 10 March, respectively « These data 
indicate a cumulative annual heating exposure of approximately 3000 Celsius 
degree-days o 
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Table 2o Summary of cumulative degree-day data and freezing exposure for six winter 
seasons 1962-63 through 1967-68. 
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Table 3o Cumulative Celsius degree-days, Milwaukee, Chicago, Two Rivers 
1962-1968 o 



1962 



1963 



1964 



1965 



1966 



1967 



1968 



Milwaukee 


6 Apro 





1 Octo 


2856 


9 Deco 


3335 


10 May 


2469 


7 Apr„ 


2608 


1 Octo 


5441 


25 Nov» 


6040 


31 May 


5533 


7 Apro 


5539 


1 Oct. 


8617 


18 NoVo 


9054 


31 Maro 


8358 


7 Apr, 


8386 


1 Octo 


11240 


16 Deco 


11718 


10 Maro 


11246 


7 Apro 


11316 


1 Octo 


14143 


25 NoVo 
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17 Maro 
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7 Apro 
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1 Octo 
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9 Dec, 
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2 Maro 
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17011 
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7216 


2 Dec. 
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31 Maro 
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14726 


6 Jan, 
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24 Febo 
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7 Apro 
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1 Oct. 
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25 Nov. 
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10 Maro 
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7 Apr. 
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1 Oct. 
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23 Dec. 


22872 


2 Mar. 


22552 


7 Apro 
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5150 
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12964 


25 Nov. 
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15189 
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The pattern of sensible heat transfer from the water mass to the atmos- 
phere of the lake may be demonstrated by taking differences between the cumu- 
lative degree-rday curves for cross-lake station pairs* Figure 3 gives the 
cross-lake difference for the northern pair of stations. Two Rivers and 
Manistee p Through the period of study, the apparent loss of heat to the 
atmosphere from the water mass as the air crosses the lake from Two Rivers 
to Manistee is equivalent to approximately 150 degree-days, between 1 October 
and 3 February as determined from air temperature measurements. This indi- 
cation cannot be interpreted in terms of a heat flux until wind data are 
considered and the volume of warmed air mass is determined. After 3 Febru- 
ary, the pattern of heat transfer becomes less distinct. The cross-lake 
warming of the air mass between Milwaukee and Muskegon is shown by the degree- 
day differences in Figure 4« For the winter seasons of 1962-63, 63-64, 64-65, 
and 66-67, the cross-lake increases were very similar (250 degree-days by 
3 March), through a range of Milwaukee freezing exposures of 507 to 866 
degree-dayso For the two seasons of mild winters 1965-66 and 67-68 (with 
freezing exposures of 472 and 439 degree-days) the cross-lake warming was 
reduced to 150 cumulative degree-days over the same period. Figure 5, 
showing the cross-lake differences for Chicago and Benton Harbor, shows no 
regular pattern of cross-lake heat exchange. In this case, it is suggested 
that the data might be biased by the southern locations of the stations, 
inasmuch as southwesterly winds do not have representative over-water tra- 
jectories between Chicago and Benton Harbor. 
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RECOMMENDATIONS 

The feasibility of using early-winter water temperature data as an ice 
potential index has been demonstrated. Following a method similar to that 
of Richards (1964), it has been shown that freezing exposure data may be 
developed to provide short-term predictions of freeze-up and break-up dates. 
Ayers' climatological calculations provide a potential means for indicating 
water temperature trends five years in advance. It is now possible to begin 
operational application of the above techniques, and to implement concentra- 
ted measurements programs to provide the necessary data for precise studies 
of the annual thermal budget of the lake. 

As suggested in a previous section of this report, all Coast Guard 
cutters should be instrumented with precision thermistor bridges for the 
measurement of intake water temperatures. Water temperature measurements 
should be made hourly during the fall runs to recover the aids to navigation. 
Additional cruises should be scheduled every three weeks during the winter 
season to obtain representative deep-water temperatures from the lake. 

Routine flights should be scheduled on a two-week basis for mapping of 
the total ice cover on the lake by both trained ice observers and using 
aerial photography as necessary o This program could be implemented by the 
Coast Guard in cooperation with the U.S« Lake Survey. 

An operational program should be established for the collection and 
maintenance of a continuous record of ice cover, winter water temperature 
data, and freezing exposure data at Coast Guard or Weather Bureau headquarters 
in Cleveland o Appropriate atlases should be compiled and maintained for 
future refinement of the operational prediction techniques. 

A cooperative program should be established between the Coast Guard, 
Weather Bureau and Lake Survey to establish solar radiation instrumentation 
at Coast Guard stations around the basin of the lake, to provide data neces^ 
sary for accurate prediction of the break-up date. 

Full support should be given to a research program in connection with 
the International Field Year Great Lakes to provide two full years of con- 
centrated study of the thermal budget of a Great Lake (Lake Ontario) . Two 
years of data will be necessary to accurately interpret the relationship 
between the lake temperature structure and the various energy exchanges 
with the air mass, radiant energy exchange, evaporation and precipitation, 
and ground-water and run-off heat exchanges. A minimum of two years of full 
data will be necessary to establish the dynamics of the thermal budget of 
the lake because of the heat storage in the water mass of the lake from 
season to seasono 
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APPENDIX As LAKE MORPHOLOGY 

Winter temperature data have been reduced and analyzed from the U.S. 
Public Health Service/Federal Water Pollution Control Administration buoy 
stations established in Lake Michigan during the winter seasons of 1962-63 
and 1963-64 o The 90-minute temperature readings from each of the 96 thermo- 
graph records have been published as Special Report Noo 39, Great Lakes 
Research Division (Noble and Michaelis 1968) o The weekly average tempera- 
tures from these records have been published in the earlier report by Heap 
and Noble (1966) o Recording thermographs were placed at 24 stations in the 
basin of Lake Michigan at depths of 10, 15, 22, 30, 60, 90^ and each succes- 
sive 30-meter level to the bottom of the lake. For completeness, and for 
future use in the computation of heat storage in the lake basin during the 
winter months, the volumes of successive depth-layers have been computed, 
using the standard thermograph station depths as centers for the depth layers. 
The depth-layers so selected were 0-45 ft, 45-67o5 ft, 67.5-87.5 ft, 87o5- 
150 ft, 250-350 ft, . 350-450 ft, etCo The areas of the top and bottom sur- 
faces of the depth layers were computed in terms of square miles, and the 
depths in feeto Accordingly the volume elements given in Table A, are ex- 
pressed in units of [mile^^ ft]e 
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Table Ao Depth-layer volumes for Lake Michigan. Layer surfaces in units 
of square miles, depths in feet. Layer volumes expressed in 
units of [mile2 ft]. 



Layer Volume 

P , L. Michigan including 

•^ Green Bay & Traverse 
Bay, [mile2 ft] 

0-45 935,131.6 

45-67.5 432,671,4 

67.5-87.5 357,621.4 

87.5-150 1,002,675.8 

150-250 1,331,181.6 

250-350 988,858.3 

350-450 567,817.6 

450-550 321,861.0 

550-650 180,948.4 

650-750 101,870.6 

750-850 38,810.2 

850-950 1,694.5 



Layer Volume 

Green Bay 

[mile2 ft] 



58,175.3 

17,850.8 

9,367.9 

2,860.0 



Layer Volume 

Traverse Bay 

[mile2 ft] 



11,381.4 
4,468.2 
3,516.8 
8,895.9 
9,231.4 
4,975.3 
1,812,2 
319.9 
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APPENDIX B^ TEMPERATURE STRUCTUJIE OF LAKE MICHIGAN 

The Great Lakes Research Division has adopted the following format for 
BT data cards i 

BT4205987170 08JU630835 0620S5005B21 000148005147007106008086010076015058021052 
BT4205987170 08JU630835 O62OS50O5B22 033048054046 

These two cards are typical of those used for BT casts o The water 
temperature is read at the water surface, and at each depth where there is 
a discontinuity in the temperature curve o These two cards are coded "BT" 
to designate the data card typeo The data are from latitude 42°05.9'N, 
87'*17oO*W, for 8 June 1963 at 0835 ESTo The data are for BT cast number 
0620 from the INLAND SEASo The BT serial number is 5005B, Two data cards 
are needed to report the cast data, and the cards shown above are cards 
numbers 1 and -2 for the given cast, respectively^ The first card shows a 
surface temperature reading of 1408*^00 The following data on the cards are 
coded as (ddd TTT), where ddd is the depth to the nearest whole meter, and 
TTT is the temperature to the nearest O^l^Co The final reading on the 
second card shows a tetnperature of 4o6®C at 54 meters. The following pro- 
gram was developed to compute the average temperature for each 10-m layer 
of water depth for each one-degree square of latitude and longitude for each 
month of the yeare 
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TEMPERATURE AVERAGES PROGRAM 





C 


THIS IS A 


PROGRAM FOR TFMPFftATlJftP AUFBACFC HP PUPRV in M nCDTU 


0001 




INTEGER DEPTH. YR.SLNO.C, BAD, PilO), IXC 30, 121, 8T/«BT«/,M0NTH(12I/ 
1»JA*,»FE«,»HA«,«AP«,'HV,"JU',«JL'.«AU'.«SP«.»OC'.«NO«,*OC«/ 


0002 
0003 


1 


REAL A(30, 121, S(300), 1110), SH(30),T£MP(121 
CALL REWIND! 71 


0004 
0005 




LSLN=0 
LCDNO=0 




0006 
0007 


2 
100 


REA0(5,10C 
F0RMATU2, 


i,END=10) ID, LT,LN,MON,YR,SLNO,N,C,IPlll,Tnj, 1-1,71 — 
A2,T4,I2,T9,i2,Ti7.A2.T19,I2.T27.l3.T36.ll.ll.T39.7II3. 


0008 




iF3.on 

IFdO.NE.BT) GO TO 2 


0009 
00 




1F(N.LT.C. 
IFtC.GT.l. 


0R.<C.EQ.1.AK0.SLN0.£Q.LSLN»I 60 TO 3 
ANO.SLNO.NE.LSLN) GO TO 3 


00 I 
0012 




IF<C.NE.tLCONO«l)) GC TC 3 
LC0NC=C 


0013 
0014 


3 


60 TO 4 
WRITE(6.101» 


0015 
0016 


101 


FURMATC0',T20. 'CARDS NCT IN ORDER') 
MRIT£(6,103)SLNO,MON,YB,N,C 


0017 
0018 


103 


F0RNAT(T20,I3,5X,A2,5X,I2,SX,I1,5X,I1) 
CALL SYSTEM 


0019 
0020 


4 


IF(C.EQ.N) 
IF(N.GE.C) 


LCDNO-0 
LSLN-SLNO 


0021 
0022 


102 


MRirE(7,102) ID,LT,Lh,MCN,YR,SLNO,N,C,(PU),T(n, 1-1,71 
FORMATIT2,A2,I2,3X,I2,6X,A2,i2,6X,I3,6X.ll,ll.(T39,7U3,F4.0»n 


0023 
0024 


10 


CALL REWINCI7) 


0025 
0026 




DO 50 LON= 
DO 50 LAT= 


85, 87 
41, 45 


0027 
0028 




00 11 K=l, 
DO 11 1=1, 


12 
30 


0029 
0030 




AI1,K)=0.0 
IX(I,K»=0 




0031 
0032 


11 

9 


CONTINUE . 
00 12 K=l, 


300 


0033 
O034 


12 


S(K)>0.0 
CONTINUE 




0035 
0036 


13 


DO 14 K=l, 
SM(K)=0.0 


30 


0037 

0038 


14 


CONTINUE 
LDEP=1 




0039 
0040 


15 


INK=0 

00 16 K=l, 


8 


0041 
0042 




P(K)=0 
T(K)=0.0 




0043 
0044 


16 
18 


CONTINUE 

READ(7,102.END=39) lO.LT.LN.HON, YR,SLNO,N,C, (P( ll.TU) , 1=1, 71 


0045 
0046 




ifilt.ne.lat.or.ln.ne.lcn.or.ioIneIbt) go TO 15 

DO 19 IM=1, 12 


0047 
0048 


19 


IFtHON.EQ. 
CONTINUE 


MONTH! IMIl GO TQ 20 


0049 
0050 


20 


H=IH 

00 21 L=l. 


8 


0051 
0052 




1F(T(1).EQ 
IFITID.EQ 


.0.01 GO TO 15 
.0.0) GO TO 22 


0053 








0054 
0055 


21 
22 


CONTINUE 
00 23 l»l, 


IN 



0056 
0057 



DEPTH-PI 11*1 
S(DEPTH)-TCII 



0058 
0059 


23 


CONTINUE 
K-LDEP 




0060 
0061 


24 


IN0-P(IN)«1 

KD-K 




0062 

0063 




00 25 K^K, INO 
IF(K.EQ.KO.OR.K.LE.IKK) GO TO 25 




0064 
0065 


25 


IF(S(KI.NE.0.0) GO TO 27 
CONTINUE 




0066 
0067 


27 


Q-<SCK0)-S1K)|/(K-K0) 
00 29 J-KO, K 




0068 
0069 




JM-J-KD 
S4J)-S«K0)-UM*Q) 




0070 
0071 


29 


CONTINUE 
IF(K.£Q.INO) GO TO 31 





0080 
0081 


33 


CONTINUE 

AtJ,MI-A(J,H) + SMU»/10. 


0082 
0083 


35 


IXtJ,M)-IXIJ,H)+l 
CONTINUE 


0084 
0085 


36 


IFIN-C) 36, 36, 37 

60 TO 9 


0086 
0087 


37 


L0EP=(IND/10I*10+1 
INK=INO 


0088 
0089 


39 


GO TO 15 

CALL REWIND(7) 


0090 
0091 


105 


MRIT£(6,105) 

FQRMATt*l*,T25, 'AVERAGE TEMPERATURE BY 10 H. LAYER') 


0092 
0093 


106 


URIT£(6,106) LAT, L(]N,YR 

FQRMATI'0',T23, 'LATITUDE-' ,T33, 12. T36, 'N LONGTITUDE-' ,T52,12, 


0094 




1TS5, 'W 19',T61,I2) 
WRITE(6,107» 


0095 


107 


FORMAT (•0*,T7, 'INTERVAL JAN fi.B MAR APR MAY JU 
IN JUL AUG SEP OCT NOV DEC) 


0096 
0097 




DO 45 1-1, 300, 10 
DO 43 J-1, 12. 1 


0098 
0099 




DEPTH- I- I 
LOEP-DEPTH+9 


0100 
0101 




IF(IX((W10)/10,J).EO.O) GO TO 42 
TEHP(J)-A((1+10)/10.J)/(10*IX((I+10)/10.J)) 


010^ 
0103 \ 


42 


GO TO 43 
TEMP(J)-0.0 


0104 
0105 


43 


CONTINUE 
IFII.EQ.1.CR.I.EQ.101.0B.I.EQ.201) GO TO 44 


0106 






0107 
0108 


108 


FORMATI* ',76, 13,110, •.•,T12,I3,T16*12(F4«i,3X)) 

GO TO 45 


0109 
0110 


44 
109 


URITEC6. 109) DEPTH, LO£P.(TEMP(J) .J«l, 12) 
FORMATt'0',T6,I3,T10,' • ,T12,13,Tie.l2IF4.l,3X» 1 


0111 
0112 


45 


CONTINUE 


0113 
0114 


110 


FORM AT fl*,T25.* NUMBER OF BT'S USED TO COMPUTE AVER. TEMPERATURE') 
WRITE(6,111) LAT, L0N,VR 


0115 


111 


FORMAT t'0*,T28,- 'LAT ITUDE-,T38,I2.T41. 'N LONGTITUDE- ,T57.12, 
1T60, 'M 19',T66,I2I 


0116 
0117 




WRiTEI6,107) 

DO 49 I-l, 300, 10 


0118 
0119 




OEPTH-I-l 
L0£P-DEPTH«9 


0120 
0121 




IFU.EQ.l.OR.I.EQ.lOl.OR.I.EQ.201) GO TO 47 
WRITEi6,113) DEPTH, LDEP, IIXU l*10)/10, J), J-l, 12» 


0122 
0123 


113 


FORMATI* •.T6,I3,T10,'_«,T12,I3,T18,12U4,3X)) 
GO TC 49 


0124 
0125 


114 
47 


FORMAT('0',T6.I3,TIO,'_',T12,13,T1«,12U4,3X)) 
WRITE(6,114) DEPTH. LDEP. < IXI ( l>10)/10.J). J-l, 12) 


0126 
01" 


49 
50 


CONTINUE 
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APPENDIX C: LAKE MICHIGAN THERMOCLINE 

The following program has been developed to objectively define the 
thermocline depth for individual BT casts • The program inspects each BT 
cast, determines the upper and lower boundary layer temperatures and depths 
for the thermocline, and defines the depth of maximum temperature gradient 
as the thermocline "meane" BT casts having only two reported readings, and 
casts with less than loO^C total temperature change are defined as isother- 
mal o In cases of complex temperature structures, two or more boundary 
layers (thermocline layers) may be defined. 
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THERMOCLINE PROGRAM 



"POB TRAM -J»-fr»COI<IM4.H^- 



— i:«T»»yi:3- <-^^-.^. - Mi^p.^^ 



THIS IS A LATITUDE BY LATITUDE MONTHLY MEAN THERMOCLINE PROGRAM 
IMPLICIT iNTEGER«4 tCtY) 

INTEQ€R*4 StNQ 

INTEGER*2 HONt MO. LMON. LAT 

DIMENSION LDI30I.LOT<30).LATt5I.M0N(12I.DEPi30).TEH130t.DELPt30l. 
10ELTI30)tTHEMT( 5*300) fDEPTP(5t3001 .THEMBI 5,300) «OEPBM{ 5«300) t 
2THTEMt5.30O).DEPTH«5.30O).ATHEMT15).ATHEM8tS).ATHTEMt5)»A0EPTPtSl. 



0010 
QOll 



0016 
00 W 



0018 
0019 



0020 
0021 



0022 
00?? 



0024 
0025 



0026 
0027 



0028 
0029 



3A0EPBM(5),A0EPTH15),NLI5) 
DATA LAT/41 ,42>»3,»4 t»S/t HGN/'JA'i 



DATA THEMr/1500*0.0/.DEPTP/1500*0.0/,DEPBH/1500*0.0/. 



ITHEMB/ 1500*0 .0/ ,THTEM/1500*0.0/, DEPTH/1500*0.0/ 
FORMAT (T4,I2tT9*I2,T17TA2,T19,12tT26,I4,T36,IltIl«T39«14I3l 



504 
50? 



FORMAT (T20, 
FORMAT (T20. 



•CARDS NOT IN ORDER AT •,T45, 14) 
'DOUBLE THERMOCLINE FOR BT «.T47.I4) 



FORMAT (T20. 'THt SECOND THERHXLINE AT • ,T47.F5.0t 5X,F5.1 ) 
FORMAT I1H0.T20. *BT' .T24. I4.T30. 'THERMOCLINE DATA ') 



FORMAT {T20, 
FORMAT tT20, 



•UPPER* 
'LOWER' 



T28,F5.0.5X,F5.1) 
.T28.F5.0.5X.F5.1) 



FORMAT (T20, 'MEAN'* T28fF5.0*5Xf F5.1) 

FORMAT (T20* 'BT'. T25.I4* T30* ' I S AN ISOTHERMAL 8T ') 



510 
511 



FORMAT UH1*TI5, 'TOTAL BT FOR'.T30, A2.T35, I2*T40* 

FORMAT CT15* 'AVERAGE THERMOCLINE DATA FOLLOW 



'IS'*T45,I3) 



514 
515 



FORMAT CT15, 'LATITUDE 4' *T25* II ,T28* 'NUMBER OF BT '*T42,I3) 
FORMAT {T15, 'DISCONTINUITY LAYER* UPPER! ' ,T45tF5.0*5X*F5.1) 



516 

520 



FORMAT (T35, 
FORMAT <T26* 



'LOWER', T45,F5.0,5X,F5.l) 

'THE THERMOCLINE IS'. T45.F5.0* 5X,F5.1) 



FORMAT ( '0', T25 

FORMAT (1H0*T10, 'TOTAL BT FOR 19'* T26*12tT30, 'IS't T35tI31 



NBT=0 
LSLNO=0 



LISO=0 

DO 1000 IM=ltl2 



DO 100 1= 
NLII)«0 



0034 
0035 




IF UM.EQ.l) GO TO 2 


00 36 
0037 


2 


GO TO 50 

DO 101 1=1. 30 


0038 
0039 




LD(I)=0 
LDTI 11=0 


0040 
0041 


101 


OELPU 1=0.0 
DELT(I) = 0.0 


0042 
0043 


3 


LCM=0 
J=1+7*LCN 


0044 


4 


K=7*7*LCN 


0046 
0047 




IF INC.EQ.CN) NBT=NBT*l 
IF <MON(IH>.eQ.MO» CO TO 5 



0048 


'5 


LMON=NO 
' HF (LSiLN0.NE.SLNO.*HU.NC.E0.CNI 60 TO 7 




0051 
0052 


6 


LCN«LCN+l 

IF (LCN.EQ.CN) GO TO 3 




0053 
0054 




WRITE! 6. 502) SLNO 
GO TO 2 




0055 
0056 


7 


LSLNO>=SLN0 




0057 
0058 




00 8 L'li 5 

IF UAT(U.EQ.LT» GO TO 9 





0106 
0107 


28 


GO TO 20 
THEMTIM.N)»sTEMin 




0108 
0109 


29 


DEPTPIH,N|sDEPII} 
JB-J 




0110 
0111 


21 


I-JB 

IF tOELPfli-D.EQ.O.O) GO TO 31 




0112 
0113 


30 


IF (OeLP(I+ll-D£LTn+I»/TOELM*tO.I 30. 30, 31 

IF ((0ELP(l4^1)/t0EPtK-l)-O£P{2t)l.GE..l) 60 TO 31 




0114 
Oil? 




T0ELH«T0ELN4-0ELT( I+l ) 
JB=JB*1 




0116 
0117 


31^ 


GO TO 21 
THEMB(M.N»*T£HII*1) 




0118 
0119 


32 


DEP8MIM,N1*DEPCI+1) 

I»JT- -....„.;, ....... 




0120 
0121 




DELPM-DEPIja*H-0£PUT) 
DEPTHIM.NMDEP(1I+0.5*0ELPH 




0122 




THT£NfN,N)<*TEN(I)^0.5*T0£LH 
WRITE! 6,505} SLNO 




0124 
0125 




WRITE! 6,506) 0EPTPCM,N) .THEMTI M,N> 
WR!TE{6.507) OEPBNIM.N)^ THENBIH.Nl 




0126 
01^7 




WR1T£{6.508) OEPTHIH,N)^ THT£HIM,N) 
GO TO 2 




0128 
0129 


33 


IS0«IS04^1 
LIS0=»LIS0+1 




0130 
0131 




URITEI6,509} SLNO 
NL1H)»NL!H)-1 




0132 
0133 


50 


GO TO 2 

WRITEt6.510) NONflM). YR. NN 




0134 
0135 




WRITE!6,550) HON! IN), YR,ISO 
WRITEI6.511) 




0136 
0137 




00 35 H-i; 5 

suMi-o:o 




0138 
0139 




SUM2>0.0 
SUH3=0.0 




0140 
0141 




SUM4»0.0 
SUH5«0.0 




0142 
0143 > 




SUH6*0.0 
y«NLIH) 




0144 
0145 




DO 34 N=l, Y 

IF (Y.EQ.O) 60 TO 34 




0146 
0147 




SUH1»SUM1+THENTIH,N) 
SUM2>SUM2+THEMBI M. N) 




0148 
0149 




SUM3»SUM3*THTEH!H*NJ 
SUN4-SUH4+DEPTPf M. Nl 




0150 
0151 




SUM5>cSUMS«0EPBMIN,N) 




0152 
0153 


34 


CONTINUE 
ATHEMTfHIsSUMl/N 




0194 
0i55 




ATHENS (N)<*SUN2/N 
ATHTEM(H)=SUN3/N 




0156 
0157 




A0EPTP(H)=SUM4/N 
ADEPBN(M)sSUM5/N 




0158 
0159 




ADEPTH(MJ»SUM6/N 
WR1TE(6.5121 M. NLl Ml 





0059 


8 
9 


CONTINUE 
M«L 


0061 
0062 




NLIM)«NLIH)«1 
00 10 I<1. 30 


0063 
0064 


10 


IF ILOT(n;£Q.O) GO TO 11 
CONTINUE 


0065 
0066 


11 


K-l-1 

IF {K.LE.3) GO TO 33 


0067 
0068 




00 13 1^1, K 
OEP!I}»LO(I) 


0069 
0070 


13 


TEM1I)«L0TII)/10. 
00 14 1=2. K 



0071 
0072 




OELTI I-l )-!TE M! I-l )-TEM! I i 

IF (0ELT!I-1).LT.0.0) GO TO 33 


0073 
0074 


14 


DELPII-1)=0EP(I1-0EP!I-1) 
N.NL(M) 


0075 




TO£LN>0.0 
T0ELS=0.0 


0077 
0078 




J-0 
Y»K-1 


0079 
0080 




DO 16 1=1, Y 

IF (TOELM.GE.OELTID) GO TO 1^ 


0081 




TOELS'TQELM 
JJ»J 


0083 
0084 




TOELM=OELT{I) 

J=I 


0085 
0086 


15 


GO TO 16 

IF (TDELS.6E.0ELTII}) GO TO 16 


0087 
0088 




TDELS=DELT(I} 
JJ=I 


0089 
0090 


16 


CONTINUE 


0091 
0092 




IF (JJ.EQ.J.0R.JJ.EQ.(J+1).0R.JJ.EQ.(J-11) 60 TO 19 
IF IT0ELS-0.5*TD£LM) 19. 18. 18 


0093 
0094 


18 


I-JJ 
DELPS=OELP(I> 


0095 
0096 




SDEP-4)EP( I )+0.5*DELPS 
STEM=TEMtI)-0.5*TDELS 


0097 
0098 




WRITE! 6*503) SLNO 


0099 
01 00 


19 
20 ) 


JT=J 
1=JT 


0101 
0102 




IF iI.LE.2) GO TO 28 


0103 
01 (14 


27 


IF I0ELPII-l)/!0EPtK-l)-0EPI2)).GE..l) GO TO 28 
TDELM«TDELM+DELTII-1) 
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